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We consider the scenario in which a damped three-level atom (in ladder- or V- configuration) is 
coupled to a single cavity mode whose vacuum state is amplified by dint of the dynamical Casimir 
effect. We obtain approximate analytical expressions and exact numerical results for the time- 
dependent probabilities, demonstrating that the presence of the third level can enrich the system 
dynamics and its population can serve as a witness of photon generation from vacuum. 



PACS numbers; 42.50.Pq, 32.80.-t, 42.50.Ct, 42.50.Hz 



O 

(N 



I 



>\ 
(N 

O 

(N 
O 
(N 



I. INTRODUCTION 

The physics of three-level quantum systems ( "atoms" ) 
interacting with quantized modes of electromagnetic field 
is very rich, and many special cases were studied in nu- 
merous papers (see, e.g., [ll-[T2| and references therein). 
In the majority of studies, the coefficients of the Hamil- 
tonians describing such systems were assumed time- 
independent. Time-dependent couplings were consid- 
ered, e.g., in [l3|, but under the restriction of adiabatic 
variations. Here we consider the light-matter dynamics 
when a three-level atom interacts with a single cavity 
mode, whose frequency is rapidly modulated according 
to the harmonical law ujt — + e sin(77i)] with a small 
modulation depth, \e\ <C 1. We shall use dimensionless 
variables, setting h = loq = 1. Such a situation can arise, 
in particular, if the selected mode describes the evolution 
of electromagnetic field in a cavity with vibrating walls, 
and one of the most impressive manifestations can be the 
so called Dynamical Casimir Effect (DCE), i.e., the pho- 
ton generation from the initial vacuum state induced by 
the motion of boundaries [13, [H. The simplest model 
describing this effect is based on the Hamiltonian [l^ 

He = ujtn-ixt {a^ - a^^) , Xt = {4,ujty^dujt/dt , 

where a and are the cavity annihilation and creation 
operators, and n = a^a is the photon number operator. 
If the modulation frequency rj is close to the parametric 
resonance frequency, rj = 2(1-1-2:) with \x\ <C 1, then one 
can expect an exponential growth of the mean number 
of photons inside the empty ideal cavity 17]. 

However, the situation can be very different if the field 
mode interacts with a detector while the cavity walls os- 
cillate. For example, it was shown in [l^ that no more 
than two photons can be created in the cavity if the 
field-atom coupling is much stronger than the modulation 
depth e. In view of the recent progress in experiments on 
simulating DCE [l9l - l2l| . the detailed study of different 
detection schemes becomes a timely and important task. 

Recently, various regimes of the two-level detectors 
were analyzed in |22h24| . On the other hand, three- 
level detectors can be more realistic [1^, besides, they 
have several advantages [26]. Therefore we consider 
the evolution of the single-mode cavity field interact- 
ing with a three-level atom, whose free Hamiltonian is 



Ha = EiUii + E21J22 + E^a^s (where Ei is the ith energy 
level, ]i) is the atomic energy eigenstate and CTy = ]i)(j]). 
We are interested in the cases where the atom-field in- 
teraction modifies severely the atomless DCE. Therefore 
we assume that the ]1) O ]2) transition is resonant with 
the unperturbed cavity frequency, E2 = Ei -\- 1. Fig. [T] 
depicts two atomic level structures we consider here: the 
ladder- (or S)- configuration and the V- configuration, 
where ^2 = E3 ~ E2 and fl^ = E3 — Ei are the transi- 
tion frequencies. We also define the detunings between 
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FIG. 1: (Color online) Diagram of the atomic configurations. 

the cavity unperturbed frequency and the other atomic 
transition frequency as A2 = 1 — ^1,2 and A3 = 1 — fJs. 
The set of constants g, g2 and 173 quantifies the atom-field 
dipolar coupling strengths between the energy levels as 
shown in Fig. [TJ We assume these constants to be real 
and 0(51) ^ 0(52) ^ 0{gz) (although this is not the 
most general case, such a choice describes the main phe- 
nomena in the most simple way) and much smaller than 
unity. Finally, the light-matter interaction Hamiltonian 
Hi must be included to obtain the total Hamiltonian 
H = He + Ha + Hi with the corresponding Schrodinger 
equation d\'i!) / dt = —iH\^). In this paper we focus on 
the amplification of the vacuum fiuctuations, so we con- 
sider only the zero-excitation initial state, ]^'(0)) — ]1,0). 

The questions we try to answer are: 1) how the pres- 
ence of the third level can influence the number of created 
photons and 2) how large can be occupation probabili- 
ties of different levels (this is important from the point 
of view of the detection of DCE). In the next two sec- 
tions we study analytically the unitary dynamics of the 
field and atom for the two configurations in the absence 
of any damping. The approximate analytical solutions 
are compared in Sec. IIVI with exact results obtained by 
solving numerically the master equation, which takes into 
account the damping in atomic degrees of freedom. How- 
ever, we neglect the dissipation in the field mode, assum- 
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ing that the cavity quahty factor is high enough. The 
last section contains a discussion of the main results. 



II. LADDER-CONFIGURATION 
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For the case depicted in Fig. [T^ we assume the Jaynes- 
Cunimings atom-field interaction Hi — 0(17(721 -l-(72f 32) + 
i/.c, where H.c. stands for the Hermitian conjugate and 
we neglected the counter-rotating terms. In the weak 
modulation case |e| 1 considered here, one can write 
Xt — 2qcos{r]t) with q = e{l + x)/4. Going to the interac- 
tion picture via the transformation |5'(t)) = V{t)\ip{t)) , 
with V {t) = exp[—it (ri/2) [n + (T33 — an)], and perform- 
ing the Rotating Wave Approximation (RWA), one ob- 
tains the approximate Hamiltonian governing the time 
evolution of |V'(i)) 

Hi ~ {gacr2i + 52acr32 — iqo^ + H.c ) 
+x[aii - n - 0-33) - A2cr33. 



where G„ = ng"^ -I- (n — 1) and An , i?„ , C„ are con- 
stants dependent on the initial conditions. Now we sub- 
stitute Eqs. ©-(EI) back into Eqs. (P)-©, assuming that 
An, Bn,Cn are slowly varying functions of time [221. If 
initially (at t = 0) all the coefficients pi„ are equal to 
zero, except for pio = 1, then pio = — (7\/2e^"'pi2 and 
for n > 2 
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A. Strong field-atom coupling 
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If the coupling between the field and atom is much 
stronger than the frequency modulation depth, |e| ^ \g\, 
we expand the wavefunction in the Fock basis as 

00 

m)) = 5]e™^*(e-*^pi„|l,n)+p2„|2,n) 



.|3,n)) 



and obtain the following differential equations for the 
probability amplitudes: 



Pin = -iVngP2{n-i) + yVi{n) 



(1) 



P2n = -igy/n + lpi(„+i) - ig2\/ne' "*P3(«-i) + W2(n) 



Pzn = -ig2\n + Te 
where for 1 = 1,2,3 



P2(„+l) + W3(n) , (2) 



Wi(n) = q v/n(n - l)e-2*^*Pi(„_2) 



-v/(r^ + l)(n + 2)e2-*pi(„+2) 



For the initial state 1 1 , 0) only the states {|l,2n), |2,2n — 
1), |3, 2n — 2)} become populated throughout the unitary 
evolution (n > 1), so the analytical analysis is simplified. 

In the resonant case, A2 = 0, the solution of Eqs. ([T])- 
© for (7 = (i.e., in the stationary cavity) is pio{t) = 
const and for n > 2 
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(A„e-^«"* - i?„e'°"* - Cn) 



(3) 



From these equations we see that by adjusting the reso- 
nance shift X [contained in the functions yVi(n)], one can 
reduce the arguments of some exponentials to zero val- 
ues, whilst other exponentials will oscillate with large 
arguments. Thus we can perform the RWA and ob- 
tain a smaller set of essential resonantly coupled differ- 
ential equations. For example, simple solutions arise for 
2x = ±G2. Then only four probabilities can be signifi- 
cantly different from zero: 



bioP ^ cos\vt), |pi2p^(^^j sin^ivt), (8) 
^ \^^^\yt), |p3o|2^i('^ysin2(^t), (9) 



where v = ^/2qg/G2 and G2 = \/2g'^ + g^- AU other 
probabilities contain extra factors of the order of (e/g)'^, 
so they can be neglected in this approximation. We see 
that at most two photons can be generated with a sig- 
nificant probability and the third level becomes partially 
populated. (If 92 = 0, then the results coincide with that 
obtained in 18].) 

A similar effect of an indirect interaction between dif- 
ferent energy levels was discovered in [23| , where the cou- 
pling constants depended on the time-dependent cavity 
length L{t) as g ^ [L(t)]~^/^, while the cavity frequency 
was supposed to be constant. One more analog is the ef- 
fect of " atomic shaking" in cavities with moving bound- 
aries, studied in [28i]. Recently, an analog of DCE in 
three-level systems with time-dependent Rabi frequen- 
cies was considered in (29l |. 

Another interesting effect is the appearance of a new 
resonance for x = (it is forbidden for the two-level 
resonant atom), due to the presence of the coefficient G2 
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in Eq. (jS]). To see this we keep only the slowly oscillating 
terms, obtaining the set of equations 



G2P10 = 2gg C2 

— iG4t ^ 



A^e-'^^* + ^26*^^* = -V6(? (^46 



52 



V2qG2P32 



and so on. We see that Pi(„>4) and P3(„>2) are cou- 
pled resonantly to C2, while C2 is coupled resonantly 
to pio, so the states |1, 2n) and |3, 2n) (n = 0, 1, 2, . . .) 
become populated throughout the evolution and many 
photons can be generated. On the other hand, the co- 
efhcients A2 and B2 are approximately uncoupled from 
Pio and from ^n>2,-Bn>2 due to the rapidly oscillating 
imaginary exponentials, so they remain approximately 
zero for all times. Therefore, the atomic state |2) is 
only virtually populated, and making the approximation 
P2n = for all times, one gets the constraint Pi(„+2) — 
— (52 /ff ) \/{n + l)/{n + 2)p3n , as we confirmed by numer- 
ical simulations. Notice that for 52 — (i.e., without the 
third level) one has C2 — therefore C2{t) = and pho- 
tons cannot be generated 18, 22]. 

In the dispersive case, IA2I » \g2\, one can write pao — 
(52/^2) e~*^^*P2i for a stationary cavity {q = 0). After 
repeating the steps above we find that for 2a; = (^2/2 ± J 
(where 62 = 92.1 A2 is the dispersive shift due to the third 
level and J = Y^(5|/4 + 2g^) the only nonzero probabili- 
ties are 



IpioP =i cos^(qi/zp<), |pi2p ~ -J^|:Sin^(qzy^t), 



(10) 



where v± = a/1 ± (52/(2 J). Thus, the rate of photon gen- 
eration and the occupation probabilities are influenced by 
the level |3), although it remains effectively unpopulated 
for all times. 



B. Weak field-atom coupling 

For a weak field-atom coupling, |e| ^ many pho- 
tons can be generated under the resonance condition 
X = 0. To calculate the accompanying atomic dynamics 
in this case, it is convenient to use the effective Hamilto- 
nian approach instead of solvin g th e differential equations 
for the probability amplitudes [2j| . We rewrite the time- 
dependent state |'(/'(0) corresponding to the Hamiltonian 
Hi as 



m)) 



1^(0)) = t/t exp i~iH,ft) t/|V(0)), (11) 



where |V'(0)) is the initial state and we introduced a 
unitary operator U to define the effective Hamiltonian 



Hef = UHiW. In the resonant regime, A2 = 0, we 
choose the transformation 24 1 



U : 



Y = (^(721 + 6o-32)+a (^(Ti2 -I- 60-23) , (12) 



where £, = 2g/e < 1 and ^ = 252/e < 1- Then to the 
second order in ^ 



Hef = iO {a^^ - a^) + iq^^2 {(Jis - cai) 



(13) 



where = q[l + ^^(ct22 - cth) + ^|(cr33 - 0-22)] is an op- 
erator with respect to the atomic basis. This effective 
Hamiltonian holds approximately for i ^ 1/g (so the 
product et can be greater than unity and several pho- 
tons can be created from vacuum). Moreover, within the 
idea of RWA we can neglect the last term in p^ . as it is 
nondiagonal in the atomic energy basis. Then the time- 
dependent wave function reads \^{t)) = C/tA(,C/|V'(0)), 
where A^, = exp[ut(a^^ — a^)] is the squeezing operator 
with nonzero matrix elements 1301 



A(«) ^ (2n|A„|0) ^ C-V2(5jc.)"^^ , 

where \k) denotes the kth Fock state and Cy = cosh {2vt), 
Sv = sinh (2vt). Using the property A^aAj = CyU — S^a^ 
one can show that for the initial state |V'(0)) = |1, 0) the 
probability amplitudes read to second order in ^ 



(1, 2nm)) = [1 - + l)]A^;;^ + e^(2n + 1)C,-/A^' 



(2,2n + IIV'(O) = i^V2^ 



e - el), while 



where Oi = q{l - C^), 6*2 = q{l 
{3,2n\ip{t)) oc ^^2, so it is much smaller and we do 
not write it here explicitly (other probability amplitudes 
are exactly zero due to the assumed initial state |1,0)). 
Therefore, in the resonant regime the third level is not 
populated within the time scale et ^ 1, but the photon 
statistics is nevertheless slightly modified due to the pres- 
ence of ^2 in 02- We checked these expressions by solving 
numerically the Schrodinger equation and found a good 
agreement. 



III. V-CONFIGURATION 

Now we repeat the steps of the previous section for 
the atom-field interaction Hamiltonian Hj = a{ga2i + 
530'3i) -|- H.c, that describes the atomic configuration 
depicted in Fig. [T|d. Using the transformation |5'(t)) = 
Vm{t)) with V{t) = exp [-ii (77/2) (n-f (733 + (722)] 
we obtain 

H2 ~ {ga<72i + g^aa^i - iqa^ + H.c.) 

-x[n + (722 + cr33) - A3Cr33. 
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For the initial state |1, 0) only the states {|1, 2n), |2, 2n — 
1), |1, 2n — 1)} become populated during the unitary evo- 
lution (n > 1). Writing the wavefunction as 

IV;) ^ ^e""*(pi„|l,n)+e">2„|2,n) 

the differential equations for the probability amplitudes 
read 

Pin = -igVnp2(n-i) - ig3Vne'*^''p3{n-i) +yVi{n) 

Pin = -igVn+ lj3i(„+i) +W2{n) 

Psn = -153 + le^^*^^Pi{n+i} + yVsin) 

and their solution in the resonant regime, A3 = 0, and 
for q — are 



Plr 



(14) 



Pl(ri-\) — -pr~ i^n^ 



B^e'^r^t „ ^ (15) 

53 



where G„ = -y^ n((7^ + gl). Therefore for weak modulation 
{\e\ <^ \g\) the resonances occur only for 2x — ±Gn, 
resulting in the probabilities 



|pio|2 ~ cos2(g<), |pi2p ^ 2 ^^=^'(9^)' (16) 



\P2l\' 



^ sm'iqt), |p3iP^(g) sin^(gt). 



(17) 

Moreover, since pi2 does not contain the coefficient C2 in 
Eq. (fT4|). the a; = resonance does not appear for the V- 
configuration. In the dispersive regime, IA3I ^ l^sj, one 

can write P31 ^ (v/253/A3)e"'*^^Pi2, where ^3 = ffl/As, 
and we find that the resonances occur for 2x = S3 ± J 
with resulting probabilities 



1 



IpioI ~ cos {qv±t), IP12I ~ sin (gft/±t) 
\P2i\ ^ ij/^sin^(gj/±i). 



(18) 



where J = -^/fJl + and i'± = ^1 ±63/ J. Notice that 
the expressions (fT8|) are slightly different from the corre- 
sponding expressions ([TUl) for the ladder-configuration. 

For strong modulation, \e\ 3> we perform the trans- 
formation 

U = e'^, r = (^(721 + e3CT3i)+a (Ccri2 + 6^113) , (19) 



where ^ = 2g/e 1 and ^3 = 2g3/e <C 1, to obtain the 
effective Hamiltonian in the resonant regime (valid for 

t « g-') 



(20) 



where fi' = g [l-f (CT22 - en 1 ) + ^| (0-33 - cri 1 )] . Neglecting 
non-diagonal atomic terms, we obtain for the initial state 
1 1,0) the following nonvanishing time dependent ampli- 
tudes: 



(l,2n|^(t)) 



[l-{e + e3){n + l)] A^^) 

+ (2n+l)[e^C,->(^)+^|C,:>^J] 



{2,2n+l\m) = iiV2^[C.^A^'^ 



A 



(")i 



-iA(«)_Ai"^l 



{3,2n+l\tP{t)) = iC3V2n:+l[Cg^A), 
where 

Oi = q{i ~e~ ^3), ^2 = qii + e), 03 = g(i + ef)- 

As expected for the resonant regime, the third level 
is populated in this case, and the photon statistics is 
severely modified as compared to the scenario of reso- 
nant two-level atom. 



IV. A COMPARISON WITH NUMERICAL 
CALCULATIONS 

We verified that our analytical results are in good 
agreement with exact numerical calculations. Here we 
demonstrate a few illustrations. To be closer to realistic 
experimental conditions, we took into account a possibil- 
ity of damping in the atomic degrees of freedom, using 
the Lindblad-type Markovian master equation for the to- 
tal statistical operator p of the atom-field system 



-i [H, p] + XD [ai2] p + Cp, 



(21) 



where Cp = X2D[a23]p for the ladder-configuration and 
— A3S)[(Ti3]p for V-configuration. The Lindblad ker- 
nel is 

D[0]p = {20pO^ - O^Op - pO^O)/2, 

A2 is the damping rate for the transition |3) — > |2), 
whereas A3 is that for the transition |3) |1). Although 
we do not consider specific atoms, we use convenient val- 
ues of g and A coefficients within the validity of our ap- 
proach and close to experimental reality. 

In Fig. [5] we show numerical results, obtained by solv- 
ing numerically the complete master equation, in regimes 
where analytical solutions are not available, for the fol- 
lowing values of dimensionless parameters: q = 3x 10~^, 
.92 =4x10-2,6 = 10-3, A = 5x10--* and A2 ^\{g2/gf. 
Such a choice agrees with a concrete example of the trans- 
mon multi-level qubit considered in [lH (assuming that 
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the noise couples to the atom via the dipolar interaction 

El). 

Nonvanishing probabilities pi„ (1 = 1,2,3) for the reso- 
nant case (A2 = 0) with the resonance shift 2x — G2 are 
shown in Fig. [2^, where we see that the third level is pop- 
ulated and at most two photons are created. The disper- 
sive case (A2 = 1032) with resonance shift 2a; = J2/2 -I- J 
is considered in Fig. ^p: as expected, the third level 
is not populated and P21 is slightly bigger than pi2, in 
accordance with the predictions given in Eq. pUj) . In 
Fig. [2t we analyze the resonance shift x ~ for the 
case of resonant third level in the absence of damping 
(A = A2 = 0). This resonance shift is absent for the res- 
onant two-level atom setup, and its dynamical behavior is 
currently not known analytically due to coupling of many 
differential equations. Still, from numerical data we see 
that for et > \ photons are generated exponentially, al- 
though with much smaller rate than for the empty-cavity 
case ((jio) denotes the mean photon number in the atom- 
less case), and the third level is significantly populated, 
while the second level remains unpopulated for all times, 
as predicted from our qualitative arguments. Thus the 
resonant ladder-configuration could be used for verifying 
photon generation from vacuum, since many photons are 
generated and the third level population can serve as the 
witness of the DCE. 
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FIG. 2: (Color online) Atom-field dynamics for the ladder- 
configuration and different resonance shifts, a) 2x = G2- b) 
2x — 82/2 + J. In (a) and (b) for et — 10 the curves are 
(from above): pio, P20 (coincident with pii), P21, P12 and pso 
[absent in (b)]. c) Behavior of average values for x = in the 
strong modulation regime. Q is the Mandel factor of the field 
mode. 

In Fig. [3] we illustrate the exact dynamics for the V- 
configuration. Fig. [3^ shows the behavior of probabilities 
in the resonant regime (A3 — 0) and weak modulation 
(e = 10^'^) for 2x — G2, while in Fig. [5]d we consider 
the dispersive regime (A3 — —12173) and the resonance 
shift 2x = 5^ + J [other parameters are: g = 3 x 10~^, 
53 = 4x10-2, A = 5x10-"*, A3 = \{g:i/gf]. In both cases 
at most two photons are created as predicted analytically 



in the absence of damping, and in the dispersive case p2i 
lies slightly above pi2 in accordance with Eq. (|18l) . In 
Fig we show the photon distribution in the absence 
of damping (obtained by tracing out the atomic degrees 
of freedom) for a; = in the strong modulation regime 
for parameters e — 10^'^, g = 5 x 10"'*, 53 = 8 x 10-"* 
and A3 = —4(73 (so the third level is neither in resonant 
nor in dispersive regime). For comparison we show the 
photon number distribution in the absence of the third 
level (^3 = 0) to emphasize that the photon statistics is 
substantially modified by the third level. For example, 
by using post-selection methods based on detecting the 
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FIG. 3: (Color online) Atom-field dynamics for the V- 
configuration and different resonance shifts, a) 2x = G2. 
For et = 10 the curves are (from above): pio, pii, P12, Pso, 
P31, P20, P21; b) 2x = 53 + J. For et = 10 the curves are 
(from above): pio, pii (coincident with P20), P21, P12. c) 
Photon number distribution for x = Q and et = 3.5 in the 
strong modulation regime and A3 = —'igz (without damping, 
A = A3 = 0), compared to the two- level atom case ((73 = 0). 



V. CONCLUSIONS 

Here we obtained simple analytical expressions for the 
atom-field dynamics when a three-level atom (see Fig. [T]) 
interacts with a single cavity field mode whose vacuum 
state is being amplified via the dynamical Casimir effect, 
and confirmed their validity by exact numerical compu- 
tations. This study is relevant since the actual atoms in 
cavity QED and artificial atoms in circuit QED are in- 
deed multi-levels. We found that the third level modifies 
the resonance frequencies as compared to the two-level 
case, and the dynamical behavior may be drastically dif- 
ferent from the cases of atomless cavity or a two-level 
atom. For instance, the modulation frequency equal to 
twice the unperturbed cavity frequency does lead to pho- 
ton generation from vacuum and occupation of the third 
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level in the ladder-configuration, whereas this modula- Acknowledgments 
tion frequency is forbidden in the case of two-level res- 
onant atom or V-configuration. Therefore the problem 
studied here might be useful to detect the Casimir pho- 
tons by measuring the third-level occupancy or creat- 
ing different entangled states. Hence the presence of the 

third level is not actually a drawback but an opportunity V.V.D. acknowledges the partial support of CNPq 
for observing a reach dynamical behavior. (Brazilian agency). 



[1] P. L. Knight and P. W. Milonni, Phys. Rep. 66, 21 
(1980). 

[2] H.-I. Yoo and J. H. Eberly, Phys. Rep. 118, 239 (1985). 
[3] S. M. Chumakov, V. V. Dodonov, and V. I. Man'ko, 

in Classical and Quantum Effects in Electrodynamics, 

edited by A. A. Komar, Proc. Lebedev Phys. Inst. Vol. 

176 (Nova Science, Commack, New York, 1988), p. 77. 
[4] V. Buzek, J. Mod. Opt. 37, 1033 (1990). 
[5] C. C. Gerry and J. H. Eberly Phys. Rev. A 42, 6805 

(1990). 

[6] Fam Le Kien and A. S. Shumovsky, Int. J. Mod. Phys. B 
5, 2287 (1991). 

[7] M. Alexanian and S. K. Bose, Phys. Rev. A 52, 2218 

(1995) . 

[8] Ying Wu and Xiaoxue Yang, Phys. Rev. A 56, 2443 
(1997). 

[9] Jing-Bo Xu and Xu-Bo Zou, Phys. Rev. A 60, 4743 
(1999). 

[10] A. B. Klimov and L. L. Sanchez-Soto, Phys. Rev. A 61, 
063802 (2000). 

[11] M. Fleischhauer and M. D. Lukin, Phys. Rev. A 65, 
022314 (2002). 

[12] A. Messina, S. Maniscalco, and A. Napoli, J. Mod. Opt. 

50, 1 (2003). 
[13] M. Janowicz, Phys. Rep. 375, 327 (2003). 
[14] V. V. Dodonov, Phys. Scr. 82, 038105 (2010). 
[15] D. A. R. Dalvit, P. A. Maia Neto, and F. D. MazziteUi, 

in Casimir Physics, edited by D. Dalvit, P. Milonni, D. 

Roberts, and F. da Rosa, Lecture Notes in Physics Vol. 

834 (Springer, Berlin, 2011), p. 419. 
[16] C. K. Law, Phys. Rev. A 49, 433 (1994). 
[17] V. V. Dodonov and A. B. Klimov, Phys. Rev. A 53, 2664 

(1996) ; G. Plunien, R. Schiitzhold, and G. Soff, Phys. 
Rev. Lett. 84, 1882 (2000); M. Crocce, D. A. R. Dalvit, 
and F. D. MazziteUi, Phys. Rev. A 64, 013808 (2001). 

[18] V. V. Dodonov, Phys. Lett. A 207, 126 (1995). 



[19] C. M. Wilson, G. Johansson, A. Pourkabirian, M. 

Simoen, J. R. Johansson, T. Duty, F. Nori, and P. Dels- 

ing, Nature (London) 479, 376 (2011). 
[20] P. D. Nation, J. R. Johansson, M. P. Blencowe, and F. 

Nori, Rev. Mod. Phys. 84, 1 (2012). 
[21] A. Agnesi, C. Braggio, G. Carugno, F. Delia Valle, G. 

Galeazzi, G. Messineo, F. Pirzio, G. Reali, and G. Ruoso, 

Rev. Sci. Instr. 82, 115107 (2011). 
[22] A. V. Dodonov, R. Lo Nardo, R. Migliore, A. Messina, 

and V. V. Dodonov, J. Phys. B 44, 225502 (2011). 
[23] A. V. Dodonov and V. V. Dodonov, Phys. Lett. A 375, 

4261 (2011). 

[24] A. V. Dodonov and V. V. Dodonov, Phys. Rev. A 85, 
015805 (2012). 

[25] B. Peropadre, G. Romero, G. Johansson, C. M. Wilson, 
E. Solano, and J. J. Garcia-RipoU, Phys. Rev. A 84, 
063834 (2011). 

[26] D. Pinotsi and A. Imamoglu, Phys. Rev. Lett. 100, 
093603 (2008). 

[27] T. Taneichi and T. Kobayashi, J. Phys. Soc. Jpn. 67, 
1594 (1998). 

[28] A. M. Fedotov, N. B. Narozhny, and Yu. E. Lozovik, 
Phys. Lett. A 274, 213 (2000); N. B. Narozhny, A. M. 
Fedotov, and Yu. E. Lozovik, Phys. Rev. A 64, 053807 
(2001). 

[29] I. Carusotto, M. Antezza, F. Bariani, S. De Liberato, and 
C. Ciuti, Phys. Rev. A 77, 063621 (2008). 

[30] R. R. Puri, Mathematical Methods of Quantum Optics 
(Springer, Berlin, 2001). 

[31] J. Koch, T. M. Yu, J. Gambetta, A. A. Houck, D. I. 
Schuster, J. Majer, A. Blais, M. H. Devoret, S. M. Girvin, 
and R. J. Schoelkopf, Phys. Rev. A 76, 042319 (2007). 

[32] M. Boissonneault, A. C. Doherty, F. R. Ong, P. 
Bertet, D. Vion, D. Esteve, and A. Blais, Preprint 
larXiv:1111.0203l 



